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Two-phase  flow  analysis 
Bolling  heat  transfer 


AM  T  A  act  Cmihm  «  «*»<>»*  'win  Ml  4r>MV  It  Hw«  •*.' 

The  power  demands  of  future  spacecraft  are  expected  to  grow  considerably  In  the  near  fu¬ 
ture.  and  this  will  require  more  efficient  thermal  transport  techniques.  Due  to  high 
heat  transfer  coefficient,  two-phase  9es-11qu1d  loops  are  expected  to  be  used  for  removal 
and  transport  of  heat  to  space  radiators.  The  design  of  such  cooling  systems  requires  a 
Knowledge  of  two-phase  fl<*  and  heat  transfer  at  reduced  and  zero  gravities.  The  objec¬ 
tive  of  this  study  was  to  develop  or  extend  the  earth  gravity  models  for  the  two-phase 
friction  multiplier,  the  void-quality  relation,  and  the  forced  convective  heat  transfer 
coefficient  to  zero  gravity  situations.  .  . 

Turbulent  mixing  length  theory  was  used  to  develop  a  model  for  the  two-phase* flow  and 
heat  transfer  parameters.  The  velocity,  density,  and  temperature  distributions  were  de¬ 
rived  by  solving  the  conservation  equations  for  a  continuous  wdlum  and  assuring  equal 
turbulent  exchanges  of  moment**,  density,  and  heat.  The  flow  was  assumed  to  be  locally 
lomogeneous, and  the  void  distribution  was  allowed  to  vary  within  the  channel.  The  resul 
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|  show  that  the  model  predicts  all  the  correct  physical  trends.  This  nodal  contains  a  free 
parameter  which  should  be  specified  from  experiments  at  zero  gravity,  simulated  zero 
gravity,  or  conditions  where  the  effect  of  gravity  can  be  neglected. 
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The  work  described  la  this  report  was  conducted  as  a  Phaas  I  Small 
tuaineas  Innovation  Research  (Slit)  program  to  begin  to  develop  better 
methods  for  calculating  two-phase  field  flow  pressure  gradieats  aed  heat 
traaafar  coefficients  that  eight  occur  in  the  weightless  eavirooeant  of 
a  spacecraft.  Nora  accurata  analytical  net hods  ere  needed  to  predict 
the  in-flight  performance  of  two-phase  fluid  heat  transfer  eyeteee  la 
comparison  to  ground  teat  results. 

The  principal  investigator  was  Dr.  Devood  Abdollahlan  with  consultation 
free  Dr.  ta lemon  Levy.  The  analytical  approach  was  based  on  earlier 
work  dome  by  the  set bore.  This  report  was  submitted  on  22  March  1585. 

Technical  direction  of  the  contract  was  provided  by  the  Plight  Dynamics 
Laboratory  of  the  Air  Pores  Mrlght  Aeronaut  leal  Laboratories  with 
Mr.  Wiliam  Raskin  serving  as  project  engineer. 

Results  of  this  effort  ere  expected  to  he  used  in  other  SRIR  contracts 
for  further  development  of  two-phase  fluid  hoot  trauafer  systems  Intend¬ 
ed  far  spacecraft  applications.  The  work  may  alao  provide  e  better 
understanding  of  basic  factors  affecting  tvo-phese  fluid  pressure  drops 
and  temperature  distributions. 
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Section  1 
MTOOOUCTION 


1.1  OACKWOUNO 

The  power  requl remen ts  for  the  future  spacecrafts  used  In  surveillance 
and  defense  missions  will  be  Increasing  as  they  9 row  In  sophistication 
and  capabilities,  military  spacecraft  are  planned  to  operate  at  power 
levels  as  high  as  a  few  hundred  kilowatts  In  the  next  20  years  (1).  All 
the  surplus  energy  penerated  In  the  spacecraft  must  be  transported  and 
ultimately  rejected  to  space.  The  only  Mode  of  heat  transfer  to  space 
Is  radiation,  and  heat  rejection  Is  accomplished  through  external 
radiators.  According  to  Reference  (2),  transport  distances  between  the 
heat  sources  and  the  radiators  vary  fron  1  to  SO  Meters,  and  transport 
capabilities  of  up  to  S  Million  watt-meters  nay  be  needed. 

The  present  thermal  control  Methods  are  based  on  transport  of  heat  via 
solid  conductors  and  internal  radiation,  heat  pipes,  and  single-phase 
liquid  and  gaseous  loops.  The  components  are  designed  for  high  tem¬ 
perature  operation  and  large  temperature  variations  when  passive  control 
techniques  are  used.  Heat  pipes  are  passive  devlres  which  can  operate 
with  small  temperature  differences  and  are  also  used  for  thermal 
transport  In  the  spacecraft.  Wien  the  solid  conductors.  Internal 
radiation,  and  heat  pipes  cannot  provide  the  necessary  transport  capabi¬ 
lity,  single-phase  fluid  loops  are  employed  In  some  spacecraft. 

The  high  operating  power  of  the  future  spacecrafts  requires  more  effec¬ 
tive  thermal  transport  techniques.  Two-phase  flow  loops  can  be  employed 
for  the  removal  and  transport  of  heat  to  space  radiators.  Mlth  this 
method,  the  operating  flow  rates  and  temperature  differences  can  be  con- 
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siderably  reduced  and  still  obtain  a  higher  heat  transfer  coefficient 
thon  a  single-phase  loop.  According  to  Reference  (2).  In  a  typical  com¬ 
parison,  the  pump  power  con  be  reduced  by  two  orders  of  magnitude  and 
the  heat  transfer  coefficient  Increased  by  two  orders  of  magnitude. 
Nuclear  power  Is  presently  being  considered  for  multi -function  com¬ 
munication  platforms,  radar  surveillance,  and  space  stations.  Two-phase 
loops  may  be  necessary  for  thermal  transport  In  the  space  nuclear  power 
systems  and  other  high  power  devices. 

Owe  to  application  In  nuclear  and  chemical  Industries,  research  In  two- 
phase  flow  has  been  In  progress  for  over  four  decades.  Even  though  It 
Is  a  very  complicated  area,  the  knowledge  In  two-phase  flow  and  heat 
transfer  has  been  greatly  advanced  In  the  last  two  decades.  Almost  all 
the  work  In  this  field  has  been  for  earth  gravity  application  and  very 
little  has  been  done  for  reduced  gravity.  Lack  of  knowledge  of  two- 
phase  behavior  at  reduced  gravities  has  prevented  Its  use  In  space 
applications.  However,  there  has  been  considerable  research  on  bubble 
growth  and  pool  boiling  which  can  be  used  for  two-phase  flew  studies. 

Experimental  and  analytical  studies  on  two  phase  flew  and  heat  transfer 
have  been  Initiated  by  MIA  and  the  Air  Force  Office  of  Scientific 
Research.  Data  and  conclusions  from  these  studies  are  needed  for  design 
of  two-phase  heat  transfer  leaps  for  space  applications. 

1.2  OBJECTIVES  AM  SMFE  V  EFFORT 

The  present  project  has  been  funded  by  Air  Force  NH0tt  Aeronautical 
Laboratories  as  a  Fhese  I  SOtR  program  to  study  two-phase  flow  and  heat 
transfer  under  aero  gravity.  The  overall  adjective  of  this  study  was  to 
develop  or  extend  the  earth  gravity  correlations  for  twe-phase  friction 
multiplier,  veld-duality  relation,  and  the  forced  convective  two-phase 
heat  transfer  coefficient  to  sere  gravity.  These  are  the  major  engi- 
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nee ring  parameters  netted  for  design  of  any  two-phase  loop.  Generally, 
the  *01 d -quality  relation  Is  needed  to  reduce  the  number  of  dependent 
variables.  The  quality  can  be  date  mined  f  row  heat  Input,  and  the 
overall  pressure  drop  can  be  calculated  by  knowing  the  void-quality 
relation.  The  twe-phase  friction  pressure  drop  Is  usually  expressed  In 
term  of  the  single-phase  pressure  drop  for  the  total  flow  considered  as 
liquid  and  a  twe-phase  Multiplier.  A  relation  for  the  two-phase 
multiplier  at  reduced  gravity  Is  needed  for  the  design  of  the  two-phase 
loops.  Aqy  device  which  deploys  the  high  heat  flux  characteristic  of 
Multiphase  flew  for  heat  transfer  will  probably  be  operating  under 
nucleate  or  forced  convective  boiling  regimes.  However,  other  flow  and 
heat  transfer  reglees  are  also  possible  for  such  equlpaents  as  conden¬ 
sers  and  boat  pipes.  Evaluation  of  the  forced  convective  heat  transfer 
coefficient  Is  essential  for  the  design  and  operation  of  two-phase  loops 
at  reduced  gravities. 

The  present  study  also  consisted  of  two  Initial  tasks  which  Included  a 
review  of  the  literature  of  earth  gravity  two-phase  flow  Models  and 
reduced  gravity  bubble  growth  Mochanlsa  and  pool  boiling.  The  analysis 
Is  limited  to  bubbly  or  slug  flow  regimes  with  no  nucleatlon  at  the 
wall.  The  heat  transfer  coefficient  Is  therefore  applicable  to  two- 
component  flows,  bubbly  flow  with  small  amount  of  vapor  at  the  wall,  or 
to  the  forced  convective  heat  transfer  regime. 
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•R1EF  REVIEW  OF  PREVIOUS  EFFORTS 

As  mentioned  In  the  previous  section,  a  review  of  the  literature  on  two- 
phase  flow  and  heat  transfer  at  earth  and  reduced  gravity  was  wade.  Our 
survey  showed  that  very  little  work  has  been  done  on  understanding  and 
modeling  two-phase  flow  at  reduced  gravities.  However,  the  areas  of 
pool  boiling  and  critical  heat  flux  have  been  addressed  considerably 
■ore.  An  excellent  review  of  the  literature  on  the  effects  of  reduced 
gravity  on  heat  transfer  for  studies  published  by  1966  Is  given  In 
Reference  (3).  Since  then,  there  has  been  wore  work  on  pool  boiling, 
bubble  growth  mechanism,  and  critical  heat  flux  but  still  very  little 
work  on  two-phase  flow.  One  reason  for  the  lack  of  activity  In  this 
area  Is  the  empirical  nature  of  aultlphase  flow  research  and  the  dif¬ 
ficulty  In  performing  such  experlaents  In  reduced  gravity.  Generally, 
the  reduced  gravity  tests  are  performed  In  laboratory  drop  towers, 
airplane  trajectory,  high  altitude  drop  tests,  and  tests  performed  In 
rockets  and  satellites.  However,  other  methods  to  simulate  reduced  gra¬ 
vity  conditions  have  been  proposed  (4),  (S),  (6),  (7),  (I).  Reference 
(3)  has  reviewed  some  flow  conditions  at  lg  where  the  effect  of  gravity 
Is  expected  to  be  negligible.  Flow  velocities  above  which  the  critical 
heat  flux  becomes  Independent  of  upward  (lg)  or  downward  (-lg)  flow 
direction  have  been  Identified. 

Results  of  an  experimental  study  (airplane  trajectory)  on  pressure  and 
temperature  changes  In  forced  convective  boiling  at  zero  gravity  have 
been  reported  In  Ref.  (9).  It  has  been  shown  that  the  system  pressure 
Increases  and  boiling  oscillations  damp  out  In  zero  gravity.  An  experi¬ 
mental  and  analytical  work  on  the  two-phase  flow  regimes  has  been 
carried  out  (10),  where  a  Froude  number  criterion  Is  Introduced  In  the 


Baker's  flow  up.  It  has  boon  concluded  that  the  pressure  drop 
Increases  significantly  as  the  gravity  decreases.  Low  heat  flux  convec¬ 
tive  boiling  at  zero  gravity  using  a  drop  tower  facility  was  studied  in 
Ref.  (11)  and.  mainly,  the  bubble  diameter  was  determined.  Forced  con¬ 
vective  maximum  heat  flux  under  reduced  gravity  was  studied  In  Ref. 
(12).  where  small  diameter  heaters  were  used  at  earth  gravity  to  simu¬ 
late  large  heaters  at  reduced  g.  Two  correlations  were  suggested  to  be 
used  under  high  and  low  flow  rates.  The  effect  of  gravity  on  the  criti¬ 
cal  heat  flux  was  studied  In  Ref.  (13).  where  tests  were  performed  in  a 
vertical  test  section  with  liquid  nitrogen  flowing  In  upward  and  down¬ 
ward  directions.  It  was  found  that  buoyancy  effect  on  the  critical  heat 
flux  decreased  with  Increasing  Inlet  velocity,  pressure,  and  subcooling. 
Other  two-phase  flow  studies  have  reviewed  the  concept  of  two-phase  loop 
application  for  spacecraft  (14)  and  (IB). 

As  mentioned  earlier,  compared  to  two-phase  flow,  considerably  more  work 
has  been  done  on  pool  boiling.  By  reviewing  analytical  and  experimental 
efforts.  Ref.  (3)  has  concluded  that  the  pool  boiling  heat  flux  Is  not 
Influenced  by  the  gravity,  while  the  critical  and  minimum  heat  fluxes 
were  found  to  decrease  with  decreasing  g.  However.  Ref.  (16)  has  shown 
that  for  liquid  hydrogen  pool  boiling,  heat  flux  Is  a  direct  function  of 
gravity.  Experiments  carried  out  by  using  a  magnetic  field  to  counter¬ 
act  the  gravitational  force,  (4)  and  (S),  have  arrived  at  the  same 
conclusion  as  Ref.  (3)  on  the  effect  of  gravity  on  nucleate  boiling  heat 
transfer.  Mechanism  of  bubble  growth  and  departure  has  been  the  subject 
of  many  Investigations  (17),  (18)  and  (6).  The  findings  are  reviewed  In 
Ref.  (2)  and  will  not  be  discussed  here. 
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Two-phase  flows  obey  the  seat  lews  of  fluid  aechanlcs  es  the  single* 
phase  flow.  There  are  basically  two  aethods  of  analysis:  differential 
and  Integral.  In  the  differential  aethod  of  analysis,  the  desired  solu¬ 
tion  Is  obtained  by  solving  the  differential  equations  which  aatheaatl- 
cally  describe  the  flow.  In  the  Integral  analysis,  the  fora  of  the 
solution  Is  asswaed  rather  than  calculated,  and  the  unknown  paraaeters 
are  date  rained  to  satisfy  the  boundary  conditions.  Integral  aethod  has 
been  used  aore  frequently  for  aodellng  two-phase  velocity  and  con¬ 
centration  profiles,  lankoff  (19),  Zuber  A  Findlay  (20). 


To  obtain  a  aathaaatlcal  nodal  for  the  flow,  two  different  approaches 
can  be  adopted.  The  two-phase  alxture  can  be  considered  as  a  single 
fluid  with  Its  own  properties,  or  the  two  phases  can  be  considered 
separately  and  conservation  equations  written  for  each  phase.  In  the 
latter  case,  the  constitutive  equations  which  describe  the  transfer  of 
aass,  nonentuo,  and  energy  across  the  Interface  oust  be  specified.  In 
the  present  analysis,  the  differential  aethod,  with  the  two  phases  con¬ 
sidered  as  a  alxture,  will  be  used.  This  approach  was  used  by  Levy  (21) 
for  evaluating  the  two-phase  density  distribution  and  pressure  drop. 
The  saae  approach  will  be  used  here,  but  the  equations  will  be  aodlfled 
for  zero  gravity  conditions  and  also  the  aethod  will  be  extended  to 
solve  for  two-phase  heat  transfer. 


If  we  consider  a  bubbly  alxture  In  zero  gravity  flowing  In  a  pipe, 
because  the  ratio  of  surface  area  to  voluae  of  the  bubbles  Is  large,  the 
drag  force  would  be  dominant,  and  the  bubbles  will  be  traveling  with 
nearly  the  saae  velocity  as  the  liquid.  Due  to  the  existence  of  trans¬ 
verse  forces,  the  bubbles  will  travel  very  close  to  but  not  exactly 
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Hong  th«  streamlines,  and  some  relative  velocity  will  exist. 
Therefore,  the  two  phases  can  be  assumed  to  flow  In  a  local  homogeneous 
state.  This  means  that  the  velocity  of  liquid  and  vapor  phases  are 
equal  at  every  radial  position,  but  there  may  be  a  concentration  gra- 
dient  In  the  pipe  cross  section. 

A  force  which  Is  normally  neglected  In  fluid  flow  analyses  but  Is  Impor- 
tant  for  void  distribution  In  zero  gravity  Is  the  lift  force  within  the 
shear  layer.  Rotation  of  a  particle  due  to  the  presence  of  a  velocity 
gradient  In  the  fluid  results  In  a  force  normal  to  the  flow  direction 
and  tends  to  move  the  particle  to  the  region  of  higher  velocity.  At 
high  relative  Reynolds  numbers,  this  phenomenon  Is  known  as  the  Magnus 
effect,  where  the  particles  experience  significant  lift  In  the  absence 
of  fluid  shear.  This  force  will  result  In  the  migration  of  the  larger 
bubbles  toward  the  center  of  the  tube  In  the  absence  of  gravity. 

3.1  METHOD  Of  ANALYSIS  AND  THE  bOVERNING  EQUATIONS 

In  the  present  analysis,  the  two-phase  system  Is  assumed  to  be  a  con¬ 
tinuous  medltm,  and  the  single-phase  turbulent  mixing  length  theory  Is 
used  to  solve  for  void  distribution  and  pressure  drop.  This  method  was 
developed  by  Levy  (21)  and  Is  extended  here  for  the  case  of  zero 
gravity. 

In  the  absence  of  gravity,  the  conservation  of  mass  and  momentum  for  a 
two-phase  flow  In  a  two-dimensional  channel  are  given  by: 

■J7  (pu)  ♦  jy  (pv)  «  0  (3-1) 


3-2 


where  x  Is  along  the  flow  direction  and  y  Is  In  the  transverse  direction 
■easured  fro*  the  wall« 

Assuming  that  the  turbulent  flow  consists  of  an  average  and  a  fluc¬ 
tuating  component ,  l.e., 

u  -  u  ♦  u'  v  •  7  ♦  v'  P  ■  ♦  P*  p  •  "p  ♦  p' 

and  taking  the  tlae  averages  of  the  conservation  equations  (3-1)  and 
(3-2),  will  result  In 

j£(e«)  (pv)  ^(eV).O  (3-3) 
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In  taking  the  tlae  averages,  the  following  averaging  rules  were  used: 


F  •  G  •  F  •  G 


F  ♦  G  ■  F  ♦  G 


The  coapresslble  boundary  layer  aoaentua  equation  Is  obtained  by 
Multiplying  equation  (3-3)  by  u  and  adding  to  equation  (3-4).  The 
resulting  equation  is 


77  Ipu2)  4  (puv)  ■  -  4  (m  ♦  jji  (-u  p'v*  -  p  uV)  (3-5) 
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Equation  (3-5)  can  be  written  In  the  following  for* 


(3-6) 


where  t  is  the  tote)  sheer  stress  which  Is  composed  of  the  viscous'  or 
leal nar  end  apparent  or  turbulent  stresses. 


t  -  tt  ♦  tt 
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(3-7) 


In  the  above  equation  It  Is  assuaeu  that  the  two-phase  viscosity  does 
not  vary  In  the  radial  direction. 

It  should  be  noted  that  the  coapresslble  turbulent  momentum  equation  In 
the  Majority  of  fluid  Mechanics  literature,  Schllchtlng  (22),  White 

(23) ,  Is  given  In  the  fora  which  wes  originally  developed  by  Van  Driest 

(24) .  This  fora  of  equation  Is  given  below. 
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and  the  total  stress  would  be 
■vjj-(MV)  (3-9) 

These  equations  are  the  saae  as  (3-4)  and  (3-7)  with  a  different  defini¬ 
tion  for  turbulent  stress.  The  definition  of  apparent  turbulent  stress, 
eq.  (3-7),  Is  Made  here  by  comparing  the  turbulent  Momentum  equation, 
(3-4)  or  (3-5),  to  the  laminar  equation.  The  reason  for  the  particular 


where  pv  ■  p  v  ♦  p'v1 
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grouping  In  equation  (3*0)  Is  that  v  appears  only  In  the  font  of  pro¬ 
duct  with  p.  However,  this  Is  not  the  case  If  a  solution  for  the  den¬ 
sity  distribution  Is  desired  when  the  grouping  of  equations  (3-0)  and 
(3-7)  Is  wore  appropriate. 

It  Is  known  that  the  temporal  average  u'v*  Is  negative  and  It  Is  nor¬ 
mally  assumed  that 


u'v*  •  -C  u*  •  v1 
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we  further  assume  that 
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Using  the  mixing  length  theory  for  the  turbulent  flow  and  absorbing  the 
constants  C  and  0  Into  the  mixing  length,  the  fluctuating  components  of 
velocities  and  density  can  be  written  as 


u*  ■  y*  ■  t  ^ 
u  ay 


(3-12) 


(3-13) 


It  Is  further  assumed  here  that  the  exchange  of  momentum  and  void  Is 
equal;  therefore,  the  mixing  lengths  for  the  velocity  and  density 
distributions  will  be  equal. 
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If  the  shear  stress  Is  approximated  by  the  wall  shear  stress,  t 

w» 

equation  (3-9)  using  the  definitions  of  equations  (3-12)  and  (3-13) 


becomes 
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(3-14) 
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dy  dy  dy  dy 


Using  tut  Van  Driest  (27)  Mixing  length  Model 


*  •  Hy  (  1  •  e  f  )l 


equation  (3-14)  beams 


tw  •  H*y2  (*)  (1  -  e  *f)2  ♦  „  ^ 

dy  dy  dy 


(3-15) 


Since  we  will  be  concerned  with  the  ween  components,  the  bar  will  not  be 
used  In  the  reMainder  of  this  section.  The  velocity,  density,  and 
pressures,  unless  specified,  will  be  the  Mean  values. 

If  another  relation  between  p  and  u  Is  found,  It  can  be  solved  slwulta- 
neously  with  equation  (3-15)  to  obtain  the  density  and  velocity  distri¬ 
butions.  In  order  to  obtain  such  a  relation,  the  forces  acting  on  the 
bubbles  will  be  analyzed  next. 

3.2  ANALYSIS  OF  THE  FORCES  ACTING  ON  A  BUBBLE  IN  ZENO  GRAVITY  TWO-PHASE 
FLOW 

The  forces  acting  on  a  bubble  during  growth  and  departure  from  a  heated 
surface  have  been  analyzed  In  References  (17)  and  (2).  These  forces  are 
buoyancy,  which  will  not  exist  In  zero  gravity;  drag;  surface  tension; 
Inertia;  and  Internal  pressure.  The  bubble  will  depart  If  the  resultant 
of  these  forces  Is  away  from  the  heated  wall.  For  Moderate  gravities. 


buoyancy  dominates  and  cause*  the  bubble  to  depart.  At  seal  1  and  zero 
gravities,  departure  My  occur  for  fluids  wh ere  the  surface  tension 
and/or  Inertia  forces  result  In  a  net  positive  force  away  froa  the 
heated  surface. 


For  flow  boiling,  a  viscous  drag  force  in  the  direction  of  the  flow  will 
be  acting  on  the  bubbles,  which  will  tend  to  detach  the  bubbles  froa  the 
wall.  For  two-phase  flow  under  zero  gravity  situations.  Internal 
pressure  and  Inertia  forces  are  negligible.  Because  the  ratio  of  the 
surface  area  to  voluae  of  the  bubbles  Is  large,  the  drag  force  will  be 
doalnant  and  the  bubbles  will  be  traveling  with  alaost  the  saae  velocity 
as  the  liquid  (local  heaogeneous  flow). 


As  aentloned  earlier,  there  Is  an  additional  force  due  to  the  rotation 
of  a  particle  which  Is  noraally  neglected  under  earth  gravity  when  the 
rotation  speed  Is  not  very  large.  This  lift  force  acts  on  a  particle 
rotating  either  due  to  the  presence  of  fluid  shear  or  rotating  freely  In 
the  absence  of  fluid  shear.  At  high  relative  velocities,  a  rotating 
sphere  will  experience  a  large  lift  force  In  the  absence  of  shear.  This 
phawonnon,  called  Magnus  effect  (25),  Is  acted  on  a  rotating  particle 
at  high  translational  velocities  like  a  golf  or  tennis  ball.  At  smII 
relative  velocities,  typical  of  situations  expected  for  two-phase  bubbly 
flow,  the  lift  force  due  to  free  particle  rotation  was  theoretically 
given  by  Rublnow  I  Keller  (25)  to  be 
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where  0  Is  the  angular  velocity  which  for  a  freely  rotating  particle  Is 
given  by 


1/2 


du 


ay 
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d  Is  the  particle  diameter,  end  p  is  the  density  of  the  fluid  around  the 
particle. 

Saffman  (27)  has  shown  that  when  the  relative  velocity  between  the  par¬ 
ticle  and  the  liquid  is  small  and  the  fluid  Reynolds  number  Is  large, 
lift  force  due  to  particle  rotation  is  less  by  an  order  of  magnitude 
than  that  due  to  shear.  The  lift  force  due  to  shear  is  given  by 

rt  *  81.2  v  P  u  d  (3-17) 


where  u  Is  the  fluid  viscosity. 

In  zero  gravity,  the  bubbles  will  be  traveling  almost  along  the 
streamlines.  Assuming  that  the  models  developed  for  particle  suspension 
Is  applicable  to  bubbles,  the  appropriate  relation  for  the  lift  force  Is 
that  given  by  equation  (3-17).  When  two  bubbles  collide,  the  resulting 
bubble  will  be  larger  and,  from  equation  (3-17),  the  lift  force  will 
grow  and  the  bubble  will  move  to  a  new  equilibrium  position  towards  the 
center  of  the  tube  which  has  a  higher  velocity. 

The  main  forces  acting  on  the  bubble  In  the  transverse  direction  are  the 
surface  tension  force,  which  holds  the  bubble  together,  and  the  lift 
force.  From  the  rules  of  the  classical  fractional  analysis,  the  ratio 
of  these  main  forces  should  be  constant  In  order  to  have  similitude  be¬ 
tween  two  systems. 
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According  to  Ref.  (28),  the  critical  or  breakup  radius  of  a  bubble  In  a 
shear  field  Is  given  by 
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(3-19) 
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It  Is  assumed  here  that  the  bubbles  coalesce  tc  the  critical  size  up  on 
Moving  to  the  new  radial  position.  Substituting  (3-19)  In  (3-18),  we 
will  get 
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The  critical  We,  according  to  Ref.  (28),  varies  fro*  0.82  to  1.0  for  the 
range  of  Ug/u^  from  zero  to  Infinity.  Therefore,  It  can  be  assumed  that 


and  eq.  (3-20)  becomes 


where  c  Is  given  by 


1  .  40.6 


(3-21) 


Dividing  both  sides  of  the  above  equation  by  the  liquid  density  and  the 
numerator  and  denominator  of  the  right-hand  side  by  tw  results  In 

u  du 

o  "  *>  ■  (3-22) 
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In  order  to  meet  the  following  boundary  conditions 
at  y  •  0  u  ■  0  o*0,  •  — 

L  ly  ii 

at  y  -  R  -  0  0  •  Pc 


equation  (3-22)  becomes 
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which  Is  the  second  relation  between  density  and  velocity  needed  for 
calculating  the  void  and  velocity  distribution. 


3.3  NOM-OIMEMSIONAl  EQUATIONS  AND  NETHOO  Of  SOLUTION 


In  order  to  solve  equations  (3-15)  and  (3-23),  they  are  non- 
dimensional  1  zed  in  terms  of  the  following  variables 
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to  get 
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£sl  -  .  1  du*  p*  + _ 2 _  (3-26) 

**'  u**>*  H2  «V**£  (1  -  *  ‘4>* 

dy  r 

The  constants  H  and  F *  for  single-phase  flo»»s  have  the  following  values: 
H  •  0.4  and  F*  -  26  (3-27) 


It  Is  known  that  these  alxlng  length  constants  are  Independent  of  fluid 
properties  and  flow  rate.  It  Is  therefore  expected  that  they  will  be 
applicable  to  turbulent  two-phase  flow.  It  Is  further  assumed  that  the 
two-phase  viscosity  can  be  approximated  by  the  liquid  viscosity.  I.e., 
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There  are  a  number  of  relations  for  two-phase  viscosity  which  can  be 
Included  In  the  present  analysis  without  complicating  the  solution. 
However,  there  Is  not  such  justification  for  applicability  of  any  of 
these  viscosity  relations.  Liquid  viscosity,  which  Is  a  good  approxlma- 
tlon  for  the  adiabatic  wall  situations.  Is  therefore  used. 


Equations  (3-25)  and  (3-26)  can  be  solved  simultaneously  for  the  two 
variables  u4  and  p4  as  a  function  of  y4.  In  order  to  have  finite  solu¬ 
tions  near  the  wall,  equation  (3-26)  should  be  Modified  for  laminar 
sublayer. 


Laminar  Sublayer  length 


Equation  (3-25)  Is  valid  for  both  the  turbulent  core  and  the  laminar 
sublayer.  However,  within  the  sublayer  u4  •  y4  and 


(3-28) 


which.  If  substituted  In  equation  (3-24),  results  In 
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The  above  equation  can  be  dl fferentlated  to  get  the  following  relation 
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At  the  edge  of  the  sublayer  y *  •  i*  the  state  variables  and  their  deri¬ 


vatives  should  be  continuous,  therefore 
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Equations  (3-30)  and  (3-31)  are  solved  for  the  sublayer  length. 
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The  above  equation  has  a  real  solution  for  pf  +  <  1/9  or  p£  <  p^/9. 

Even  thou9h  the  two  non -dimensional  state  variables  u*  and  p+  can  be 
solved  as  a  function  of  y*  from  the  above  equations,  the  results  Mill 
not  have  much  practical  use  unless  rearranged  In  other  forms.  In  order 
to  do  this,  the  following  relations  are  derived. 


The  cross-sectional  average  density  Is  given  by 


Considering  that  r  •  R  -  y  and  dr  ■  -dy,  the  above  equation  can  be  non- 
dimensional  1  zed  to  get 
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where  Q  Is  a  dummy  variable. 


The  average  mass  flux  Is  defined  as 
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This  equation  can  also  be  non-dlmenslonallzed  as  follows  to  calculate  a 
new  variable,  2. 
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This  two-phase  Reynolds  nuwber  Is  defined  as 

Re  .  £!£*!  «  2  H*  2  (3-35) 
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The  two-phase  friction  ■ultlpller  Is  defined  as 


(3-35) 


where  the  denoalnator  In  the  above  equation  Is  the  frictional  pressure 
gradient  If  the  total  flow  Is  assumed  to  flow  as  liquid,  l.e., 
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where  fQ  Is  the  single  phase  liquid  friction  coefficient  given  by  the 
following  relations: 


f  ■  —  Re  <  2300 

0  Re 

f0  •  1.84  x  10‘4  Re0**44  2300  <  Re  <  4000  (3-38) 
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fQ  •  0.316  Re" 


Re  >  4000 


The  we1!  sheer  Is  given  by 


Therefore,  two-phese  friction  multiplier  becomes 
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The  void  frectlon  for  the  two-phese  flow  Is  celculeted  from 
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The  eree  evereged  flow  quellty,  Is  obtetned  from  the  following 

re  let  ion  for  locel  homogeneous  flow 

“g 

%  m  -  .  - 

*  W  GA 


WG  "  PG  *  “G  *  °G  *  u 


3-16 


Non-dimensional  1  ling  the  equation  results  In 


(3-41) 


Using  equations  (3-34)  and  (3-41),  the  flow  quality  Is  calculated  as 
follows: 


*a 


(lV)  u*  (**  -  0)  do 


P t  u4  (N4  -  0)  d0 


(3-42) 


The  cross-sectional  area  average  void  becomes 


1 


% 


He  next  define  the  following  void  and  quality  relations  to  absorb  the 
effect  of  system  pressure. 
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(lV)  u*  (R*  -  Q)  dO 
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(3-44) 


Method  of  Solution 

for  a  given  value  of  the  constant  c: 

1.  A  non-d 1 Mens  1 ona 1 1  zed  center  line  density,  pc*.  Is  assuMed. 

2.  The  non-d Inenslonal  sublayer  length,  6*,  Is  calculated  from 
equation  (3-32). 

3.  The  governing  equations  are  Integrated  using  equations  (3-28)  and 
(3-29)  for  y*  <  4*  and  equations  (3-2S)  and  (3-26)  for  y*  >  t*  to 
calculate  the  state  variables  u*  and  p*  as  a  function  of  y* . 

4.  At  every  step,  the  Integrals  (3-33),  (3-34),  and  (3-44)  are  per- 

forned  to  calculate  p  ♦  2  4 

’  *  *o  •  “a  *  *nd  xa*  fro*  equations 
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(3-35) ,  (3-39),  (3-43),  ami  (3-44) 


Therefore,  for  a  given  value  of  the  constant  c,  the  state  variables  are 
calculated  as  a  function  of  Y*  and  p£4.  In  addition,  because  for  every 
y*  and  *>c+,  there  are  specific  Re,  p^4,  *«d  ,  the  void 
fraction  and  two-phase  friction  Mltlpller  can  be  expressed  as  a  func¬ 
tion  of  x4*  ond  R«. 

A  computer  prog raw,  ZROG,  was  developed  to  perfom  the  above  calcula¬ 
tions.  This  program  uses  the  SFOOC  Integration  algorithm  which  solves  a 
system  of  ordinary  differential  equations  using  Adams  predictor- 
corrector  method.  ZROG  was  written  In  FORTRAN,  and  the  calculations 
were  performed  on  a  Compupro  (with  NP/N  operating  system)  and  on  an 
IBM- PC/XT  (with  MS-DOS  operating  system).  All  the  routines  were  written 
In  double  precision  to  reduce  the  round-off  errors.  The  program  GRAFIT 
(developed  by  Golden  Software)  was  used  to  generate  the  plots.  This 
routine  only  connects  straight  lines  between  points  and  no  curve  fitting 
Is  performed. 

3.4  RESULTS  AMO  DISCUSSION 

The  parameter  c  Introduced  In  equation  (3-21)  Is  actually  a  propor¬ 
tionality  constant  which  Is  a  combination  of  the  constant  X  and  the 
similitude  ratio  K^. 

C  .LUl.tSJ. 

*3  *3 

The  value  of  c  should  be  determined  from  experiments,  but  because  It 
affects  the  length  of  the  laminar  sublayer  given  by  equation  (3-32), 
some  order  or  magnitude  approximation  for  this  variable  can  be  made. 

Van  Drlest's  (24)  mixing  length  equation  used  here  Is  actually  a  con¬ 
tinuous  relation  which  can  provide  an  eddy  dlffuslvlty  applicable  all 


the  way  to  the  wall.  However,  equations  (3-25)  and  (3*26)  are  not  con* 
tlnuous  close  to  the  wall,  trtilch  aates  It  necessary  to  Introduce  a 
sublayer  region,  as  discussed  before,  lecaute  the  aolecular  dtffuslvity 
Is  expected  to  die  out  for  y4  -  10  to  30,  fro*  Table  3*1,  which  Is  for 
pc4  -  0.1,  a  value  of  the  order  of  1/100  sews  to  be  appropriate.  Table 
3.2  gives  the  sublayer  length  for  c  •  1/200  and  1/400  for  different 
values  of  centerline  density. 


Table  3-1 

I4  Values  for  •  0.1 


c 

1 

1/10 

1/100 

1/200 

1/400 

1/1000 

«4 

1.49 

4.71 

14.91 

21.06 

29,  tl 

47.14 

Table  3*2 

4 4  Values  for  c  •  1/200  and  1/400  *- 


'c* 

0.01 

0.02 

n 

jm 

B 

B 

B 

B 

B 

0.1 

c-1/200 

14.51 

14.91 

15.34 

15.83 

16.36 

16.97 

17.68 

18.52 

19.59 

21.06 

C-l/400 

20.52 

21.06 

21.70 

22.38 

23.14 

24.0 

25.0 

26.19 

27.7 

29.81 

Velocity  and  density  distributions  for  pc*  .  QA  #n„  for  differtnt 
values  of  c  are  shown  In  Figures  3.1  and  3.2.  The  density  profile  be- 
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coats  flitter  for  smaller  c.  Figures  3.3  to  3.6  show  u*  and  p*/pc*  for 
c  •  1/200  and  c  •  1/400.  As  discussed  In  the  previous  section,  for  the 
given  value  of  c,  two-phase  friction  Multiplier  and  void  fraction  can  be 
expressed  as  a  function  of  flow  quality  and  Reynolds  number.  Figures 
3.7  to  3.10  show  40*f#  and  a*  as  a  function  of  **  for  c  ■  1/200  and 
c  •  1/400.  Predictions  are  presented  In  this  fora  to  absorb  the  effect 
of  system  pressure  and  the  choice  of  single-phase  friction  coefficient 
relation. 

The  two-phase  friction  multiplier  using  the  single-phase  friction  factor 
selected  here,  equation  (3-38),  Is  shown  In  Figures  3.11  and  3.12.  At 
earth  gravity,  the  two-phase  friction  multiplier  Is  expected  to  decrease 
with  Increasing  mass  flux.  An  opposite  trend  Is  shown  <n  Figure  3.11 
and  3.12.  This  Is  due  to  the  fact  that,  even  though  40^f0  drops  with 
Increasing  auss  flux,  fQ  given  by  equation  (3-381  has  a  stronger  depen¬ 
dence  on  Re  which  results  In  an  Increasing  two-phase  multiplier.  For 
larger  values  of  c,  the  two-phase  multiplier  My  show  the  same  trend  as 
the  earth  gravity  case.  In  any  event,  the  true  coefficient  for  calcu¬ 
lating  the  two-phase  pressure  drop,  as  shown  In  equations  (3-36)  and 
(3-37)  Is  4C  fQ,  which  reduces  with  Increasing  mss  flow.  In  Figure 
3.11,  the  two-phase  friction  multiplier  for  Reynolds  number  of  3000  Is 
larger  than  for  Reynolds  number  of  S000.  The  reason  for  this  trend  Is 
that  the  single-phase  friction  coefficient  In  the  transition  region  be¬ 
tween  Re  •  2300  and  Re  ■  4000  goes  down  which  results  In  an  Increasing 
two-phase  friction  multiplier.  As  mentioned  earlier,  the  ^lotting 
routine  used  here  connects  the  data  points  with  straight  lines. 
Therefore,  slope  changes  such  as  the  ones  seer  on  Figures  3.2  and  3.6 
around  y4«10  could  have  beer  avoided  If  more  data  points  were  available 
or  a  curve  smoothing  technique  was  used. 

The  predictions  using  the  present  model  can  be  compared  to  earth  gravity 
data.  Unfortunately,  the  void-quality  data  In  the  Mjorlty  of  the  two- 


phase  flow  literature  Is  presented  either  as  a  function  of  thenaodynaatlc 
quality  or  it  Is  reported  for  very  low  qualities.  The  void  as  a  func¬ 
tion  of  flow  quality  can  be  obtained  fro*  the  available  literature,  but 
this  requires  more  time.  Due  to  project  Imitations,  this  coaparlson  Is 
not  perforaed  In  this  phase  of  the  program.  However,  the  two-phase 
friction  aultlpller  froa  the  present  aodel  Is  coapared  to  earth  gravity 
data.  This  Is  shown  In  figure  3.13  which  Is  reproduced  froa  Reference 
(29).  The  two-phase  friction  aultlpller  for  Re  *  50,000  and  c  ■  1/200 
and  1/400  are  shown  as  a  function  of  flow  quality.  Froa  Figures  3.7  to 
3.13,  It  can  be  concluded  that  the  paraaeter  c  does  not  have  a  signifi¬ 
cant  effect  on  the  void-quality  relation,  while  It  strongly  affects  the 
two-phase  friction  multiplier.  Figure  3.13  also  shows  that  the  two- 
phase  aultlpller  increases  with  Increasing  the  paraaeter  c  and,  for 
values  of  1/200  and  1/400,  the  two-phase  aultlpller  Is  significantly 
saaller  than  the  earth  gravity  results.  This  Is  In  contrast  to  conclu¬ 
sions  of  Ref.  (10),  where  the  pressure  drop  was  shown  to  Increase  as  the 
gravity  was  reduced.  It  Is  believed  that  the  present  model  represents 
the  physical  phenomenon  correctly,  but  the  value  of  c  used  In  the  pre¬ 
dictions  aay  be  too  small.  As  discussed  In  Section  5,  this  aodel  pre¬ 
dicts  an  Increasing  void  fraction  approaching  the  homogeneous 
distribution  and  a  decreasing  two-phase  friction  coefficient,  40  fQ,  as 
Reynolds  number  Is  Increased.  The  two-phase  multiplier  should  be  larger 
than,  but  approaching  the  homogeneous  distribution  as  the  Reynolds 
number  Is  Increased.  Figure  3.13  shows  that  the  two-phase  multiplier  Is 
smaller  than  the  homogeneous  value.  A  larger  value  of  c  (on  the  order 
of  1/20  to  1/50)  would  result  In  a  larger  two-phase  multiplier  and  a 
decreasing  trend  with  increasing  Reynolds  number.  Such  values  of  c 
correspond  to  sublayer  lengths  of  5  to  10  and  will  not  affect  the  void- 
quality  relation  significantly.  Due  to  time  limitations,  predictions 
with  larger  values  of  c  were  not  performed  In  this  phase  of  the  project 
but  are  recommended  for  the  future. 


-  - - 
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Following  the  drift  flu*  approach  (20).  the  generalized  void-quality 
Model  can  be  expressed  as 


[  x  ♦  —  (1-x)  ]  4 

pl 


ocV 


(3-45) 


where  cQ  Is  the  void  distribution  parameter  and  Is  the  drift  velo¬ 
city.  For  the  case  of  local  homogeneous  flow,  »  0  and 

«  .  - 1 -  (3-46) 

t  x  *  —  (1-x)  ) 

°L 

cQ  Is  usually  determined  fro*  experiments,  and  It  varies  between  1.0  and 
1.2  for  earth  gravity  bubbly  flow.  In  the  present  analysis,  an  attempt 
was  made  to  find  the  void  distribution  analytically  using  the  analysis 
of  the  forces  acting  on  the  bubbles.  However,  there  Is  an  unknown  fac¬ 
tor,  c,  which  should  be  determined  from  experiments.  At  very  high 
pressures,  the  two-phase  flow  at  earth  gravity  will  be  nearly  homoge¬ 
neous.  In  addition,  decreases  with  Increasing  pressure,  thereby 

reducing  the  effect  of  gravity  on  the  flow.  In  the  absence  of  zero  gra¬ 
vity  test  results,  high  pressure  vertical  two-phase  flow  tests  can  be 
used  to  estimate  the  parameter  c. 


LOG  r+ 
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Figure  3.2.  Efftct  of  the  Parameter  c  on  the  Density  Profile 
for  »  0.1 
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Figure  3.5.  Effect  of  Center  Line  Oenslty  on  Velocity  Profile 
for  C  •  1/400 
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Figure  3.7.  Two-Phase  Friction  Factor  vs.  Modified  Flow  Quality 
for  C  •  1/200 
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Figure  3.8.  Void-Quality  Profile  for  C  *  1/200 
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Figure  3.9.  Two-Phase  Friction  Factor  vs.  Modified  Flow  Quality 
for  C  •  1/400 


□  Ri-3000  ART  =5000  ORE*  10000  XRE*20000  +RE*SOOOO 

Figure  3.11.  Two-Phase  Friction  Multiplier  vs  Modified  Flow  Quality 
for  C  »  1/200 
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Section  4 

TWO-PHASE  HEAT  TRANSFER  COEFFICIENT 

As  discussed  In  Section  2,  the  two-phase  Neat  transfer  coefficient  at 
earth  gravity  depends  on  the  flow  regime  and  whether  the  bulk  liquid  Is 
subcooled  or  saturated.  The  heat  transfer  regimes  for  forced  convective 
boiling  are  shown  In  Figure  4.1.  The  boundaries  of  the  saturated  forced 
convective  boiling  are  shown  In  Figure  4.2.  The  present  analysis  is  for 
saturated  two-phase  heat  transfer  which  covers  the  nucleate  boiling 
regime  (regions  C  and  0  In  Figure  4.1)  and  the  forced  convective  regime 
(regions  E  and  F  In  Figure  4.1)  where  there  Is  no  nucleatlon  at  the 
wall.  As  shown  In  Figure  4.2,  the  forced  convective  boiling  can  occur 
at  any  flow  pattern,  but  normally  at  earth  gravity  this  heat  transfer 
regime  occurs  with  the  annular  flow. 

The  majority  of  the  existing  correlations  (30)  and  (31)  are  of  the  form 


Np  B  hfo  r 

where  h^p  and  hfQ  are  the  two-phase  and  single-phase  liquid  heat 
transfer  coefficients,  respectively,  and  x-t  Is  Hartlnelll  parameter. 
These  correlations  have  been  developed  for  the  forced  convective  region 
and  have  been  modified  for  the  application  to  saturated  nucleate  boiling 
regime.  The  Chen  (32)  correlation,  however.  Is  applicable  to  both 
regimes  and  has  been  developed  by  adding  the  contribution  due  to  both 
mechanisms  of  heat  transfer. 

The  present  analysis  should  basically  be  valid  to  any  mode  of  heat 
transfer  which  Is  controlled  by  the  temperature  difference  In  a  narrow 
region  close  to  the  wall.  However,  because  the  velocity  and  density 
distributions  developed  In  the  previous  sectlor  are  used,  the  analysis 
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Is  valid  for  no  vapor  bubble  at  the  wall.  This  will  cover  heat  transfer 
for  two-phase,  two-component  flow  and  single- component  forced  convection 
regime  with  bubbles  within  the  liquid  film  or  for  nucleate  boiling  with 
a  mall  number  of  oucleatlor  sites,  where  fluid  density  at  the  wall  can 
be  approximated  by  the  liquid  density. 

4.1  GOVERNING  EQUATIONS 

^e  conservation  of  energy  for  a  two-dimensional  channel  Is  given  by 

oc  (u  ♦  ¥  il)  -  -i  (K  il)  ♦  u  ^  ♦  up  (4-1) 

¥  >x  ay  ay  ay  dx 

where  a  Is  the  dissipation  term.  Following  the  same  approach  as  in  the 
previous  section,  it  can  be  assumed  that  for  the  two-phase  turbulent 
flow,  the  temperature  consists  of  an  average  and  a  fluctuating  com¬ 
ponent,  l.e., 

T  «  T  ♦  T* 

Using  the  averaging  rules  shown  In  Section  2  and  taking  the  time  average 
of  equation  (4-1),  we  will  get 


♦  C  0  V  — 

P  ay 


(«  £>  -  c_  eV  il  -  (c,;  v'T' ) 

ay  ay  P  ay  ay  P 


♦  U  &  ♦  U# 
dx 

If  It  Is  assumed  that  the  specific  heat  does  not  vary  with  y,  the  above 
equation  reduces  to 


c  p  u  —  ♦  c  p  v  —  ■  — i  (K  — )  -  c  p*v*  il  -  -i  (p  v’T* ) 
"  ax  P  ay  )y  ay  P  ay  ay 


(«-2) 


♦  U  ^  ♦  U# 
dx 


In  the  compressible  flow  literature,  the  tern  -  p'v‘  aT/ay  In  the 
above  equation  1$  combined  with  the  transverse  convective  tens  to  get 

(4-3) 


c.  c„  ,7il.  -i3»2£.  * 

p  ax  p  ay  ay  ax 


where  In  comparison  to  laminar  energy  equation 


q  "  \am1nar  4  ^turbulent 

q  ’  '  K  fy  *  CP  *  *'"**  («-*> 


The  above  definition  for  the  heat  flux  Is  equivalent  to  assuming 

®* v*  T  <  <  v'T*  p  (4-5) 

In  equation  (4-2)  and  defining  (4-4)  by  comparing  (4-2)  to  laminar 
energy  equation. 

In  the  superheated  region  near  the  wall.  It  1$  assumed  that 
«  *  «w 
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Applying  the  Mixing  length  theory  end  using  equal  Mixing  lengths  for  the 
MOMentiM  end  he«t  transfer,  we  Mill  get 


«w 


s  HV  c 
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(4-5) 


4.2  NON-OIHENSIONAI  EQUATIONS  AM)  NETHOO  OF  SOLUTION 

Equation  (4-5)  Is  non-dlwenslonallzed  using  the  definitions  of  equation 
(3-23)  and 

*  C«PL  pl 

T  -  (T  -  T)  (4-5) 

«w 


If  we  drop  the  bar  convention  representing  the  wean  component  of  the 
flow,  (4-5)  reduces  to 
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du  4  _1 
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dy 


(4-7) 


It  Is  known  that  the  two-phase  Prandtl  nuMber  can  be  approxlnated  by  the 
value  corresponding  to  the  liquid,  l.e., 

Pr  •  Pr^ 

Equation  (4-7)  can  be  Integrated  along  with  the  differential  equations 
for  p*  and  u4  given  by  (3-29)  and  (3-2f)  for  y4  <  l4  and  equations 
(3-24)  and  (3-25)  for  y4  >  44.  Therefore,  for  a  given  value  of  the 
constant  c,  the  state  variable  T4  will  be  calculated  as  a  function  of 
y *  and  a,* 

V  M 
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In  order  to  calculate  the  heat  transfer  coefficient,  a  mixed  mean  tem¬ 
perature  Is  defined  as  follows 


6  A  cp  Tm  " J  pu  cp  T  dA 


For  constant  c„,  we  will  have 
P 


(4-8) 
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ton-dimensional Izln?  (4-9)  and  defining  a  heat  transfer  coefficient  as 


(«-10) 


will  result  In 


T  *  ^  W*«r2  f*  **  r*  I  m++\M*  1 

T.  •  T.  *  t  *  .  •«  T  («-»)<»  1 


^  __ 

ij  I"  .V»*  <«V»  «r‘ 


(4-11) 


.•  .*  .*■  *• 


Nusselt  number  1$  defined  as 
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The  aethod  of  solution  consists  of  Integrating  (4-7)  along  with 
equations  (3-28)  and  (3-29)  for  y4  <  4*  and  equations  (3-25)  and  (3-26) 
for  y4  >  «4 .  At  every  step  the  Integral  In  equation  (4-12)  Is  performed 
and  a  Nusselt  number  Is  calculated.  All  the  solutions  are  perforated  for 
Pr  ■  1.0  In  the  present  analysis. 

The  computer  program  ZROfc,  discussed  earlier,  actually  solves  the  system 
of  three  differential  equations  for  the  state  variables  u4,  p4,  and  T4. 


4.3  RESULTS  AND  DISCUSSION 

Temperature  distributions  for  a  given  value  of  the  center  line  density 
(p£4  *  0.1)  and  different  values  of  c  are  shown  In  Figure  4.3.  The 
change  In  the  slope  In  this  figure  at  around  y4  •  10  Is  not  really  pre¬ 
dicted  by  the  model.  This  Is  due  to  logarithmic  scales,  the  few 
available  points  between  y4  •  0  and  10,  and  the  Method  of  plotting. 
However,  the  change  In  slope  and  Increase  In  temperature  at  larger 
values  of  y4  (e.g.,  y4  •  100  for  c  ■  1/400)  Is  representative  of  the 
true  phenomenon.  This  rise  In  temperature  at  larger  y4  is  due  to 
Increased  void  within  the  flow.  Figures  4.4  and  4.5  show  the  tem¬ 
perature  distribution  for  c  ■  1/20U  and  1/400  and  different  values  of 
the  center  line  density. 


The  Nusselt  number  predicted  by  equation  (4-12)  Is  not  a  function  of  Re 
4hd  P  r^  alone,  and  It  varies  with  the  center  line  density.  Similar 
solutions  for  the  single-phase  turbulent  flow  In  a  pipe  express  Nu  In 
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terns  of  Re  end  Rr  alone.  The  reason  for  this  Is  that.  In  addition  to 
the  differential  Method  as  described  for  the  two-phase  flow  above,  the 
single-phase  Models  assuee  a  velocity  profile  to  calculate  the  state 
variables  *t  the  tube  center  line. 

The  husselt  metiers  for  c  ■  1/200  and  1/400  for  different  Re  and  pc*  are 
shown  in  Tables  4-1  and  4-2.  It  can  be  seen  that  Hu  Increases  with  Re 
as  expected.  Hu  decreases  with  Increasing  center  line  density  but 
reverses  this  trend  at  certain  Reynolds  nwnbers.  The  reason  for  this 
reversal  In  the  trend  needs  to  be  Investigated. 
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Figure  4.3.  Effect  of  the  Parameter  c  on  the  Temperature  Profile 
for  oc*  •  0.1 
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Section  5 


conclusions  and  recommendations 


Turbulent  mixing  length  theory  was  used  to  develop  models  for  two-phase 
friction  multiplier,  void-quality  relation,  and  two-phase  heat  transfer 
coefficient  under  zero  or  negligible  gravity.  Differential  method  of 
approach  was  used  to  solve  the  governing  equations  for  the  conservation 
of  mass,  momentum,  and  energy.  The  mixing  length  model  of  Van  Driest 
(24)  was  used  here,  and  a  sublayer  region  was  considered  to  avoid 
discontinuity  In  the  numerical  solution.  The  non-dimensional  velocity, 
density,  and  temperature  were  used  to  solve  for  the  desired  variables. 
This  model  contains  a  free  parameter,  c,  which  Ideally  should  be  deter¬ 
mined  from  experimental  results. 

In  this  analysis,  the  flow  was  assumed  to  be  locally  homogeneous,  but 
the  void  distribution  was  allowed  to  vary  within  the  channel  cross  sec¬ 
tion.  The  driving  force  for  this  distribution  Is  a  lift  force  which 
acts  on  the  bubbles  and  drives  them  to  the  region  of  higher  velocity. 
Physically,  this  model  predicts  the  correct  trends  for  the  two-phase 
flow.  As  the  Reynolds  number  Increases: 

•  The  bubbles  will  be  more  uniformly  distributed  and  the  flow  becomes 
more  homogeneous. 

•  The  void  fraction  Increases  and  the  distribution  approaches  the 
homogeneous  distribution. 

•  The  two-phase  friction  factor,  a  2  f  .  decreases. 

D  0 

•  Nusselt  number  Increases. 

However,  It  Is  believed  that  the  values  for  c  used  In  the  predictions 
are  too  small  and  c  on  the  order  of  1/20  to  1/50  will  be  more 


appropriate.  This  will  result  In  two-phase  multipliers  larger  than 
homogeneous  which  will  approach  homogeneous  distribution  as  Reynolds 
number  Is  Increased. 

The  Nusselt  number  predicted  by  this  model  Is  a  function  of  Reynolds 
number,  Prandtl  number,  and  the  velocity  and  density  distributions  which 
are  specified  through  the  center  line  density.  The  predicted  Nusselt 
numbers  reverse  their  dependence  on  center  line  density  at  a  specific 
Reynolds  number,  as  shown  In  Tables  4.1  and  4.2.  This  Is  related  to  the 
temperature  distribution  as  shown  In  Figures  4.3  to  4.9. 

Certain  trends  predicted  by  this  model  need  to  be  Investigated  further 
to  establish  the  theoretical  bases  of  the  model.  These  are: 

1.  predictions  should  be  performed  for  larger  values  of  c  (1/20  to 
1/100)  and  the  results  compared  to  homogeneous  distribution  as  flow 
velocity  Increases. 

2.  The  assumption  of  r  •  and  q  -  should  be  evaluated, and  the 
possibility  of  Including  the  exact  relations 

T  "  Tw  d  *  R  *  S»  U  -  £) 

n  n 

should  be  Investigated. 

3.  The  void-quality  and  two-phase  multiplier  should  be  compared  to 
more  earth  gravity  data. 

4.  The  temperature  profile  and  its  dependence  on  center  density  should 
be  studied. 

5.  The  effect  of  Rrandtl  number  on  temperature  distribution  and 
Nusselt  number  should  be  Investigated. 


6.  HvtlMtion  at  the  wall  was  neglected  In  the  present  model.  The 
possibility  of  extending  the  model  to  Include  nucleatlon  needs  to 
be  Investigated. 

The  free  parameter  In  the  present  model  needs  to  be  specified  aiore 
accurately.  Obviously,  experimental  results  under  reduced  gravity  or 
conditions  which  simulate  the  flow  behavior  at  zero  gravity  are  the 
ultimate  test  of  the  model.  These  types  of  tests  should  be  considered 
In  the  future,  and  they  can  be  used  to  establish  the  parameter  c.  In 
the  absence  of  such  experiments,  high  pressure  two-phase  flow  test 
resu'ts  can  be  used  for  better  estimation  of  c.  An  alternative  approach 
should  also  be  considered  which  consists  of  selecting  a  constant 
sublayer  length  and  allowing  the  model  to  calculate  c,  velocity,  den- 
slty.  and  temperature  distributions.  A  series  of  two-phase  flow  tests 
have  been  designed  and  planned  for  the  future  space  shuttle  flights 
(36).  The  design  for  these  experiments  Is  based  on  earth  gravity 
models.  It  will  be  very  beneficial  to  review  these  designs  In  view  of 
the  predictions  by  the  present  model  and  also  study  the  possibility  of 
enhancing  these  tests  to  calculate  the  parameters  needed  for  this  model, 
which  seems  to  be  the  only  one  for  zero  gravity. 

In  addition  to  the  above  recommendations,  which  are  mainly  for 
establishing  the  present  approach,  efforts  are  needed  to  study  the  flow 
regime  transitions  at  zero  gravity  and  also  other  possible  two-phase 
flow  and  heat  transfer  regimes. 
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